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Abstract 

We propose to generalize Bekenstein model for the time variation 
of the fine structure "constant" a em to QCD strong coupling constant 
as- We find that, except for a "fine tuned" choice of the free pa- 
rameters, the extension can not be performed trivially without being 
in conflict with experimental constraints and this rules out as vari- 
ability. This is due largely to the huge numerical value of the QCD 
vacuum gluon condensate when compared to the matter density of the 
universe. 
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1 Introduction 



The time variation of fundamental constants may provide a connection be- 
tween cosmology and particle physics. Early suggestions can be traced to 
Dirac [0 long ago, but many proposals leading to time-varying constants 
were discussed far afterwards. These can be classified to either "phenomeno- 
logical" models 0, 0, J§, models providing a natural theoretical framework 
in terms of higher-dimensional theories like Kaluza-Klein || || and string 
theories |7|, |3| or models based on first principles such as the Bekenstein's 
for a em variability || . Besides, spatial variations of the fine structure within 
Beckenstein's models have also been considered |TD| 



Actually, Bekenstein model for electromagnetism is very attractive be- 
cause it is based on very general assumptions: covariance, gauge invariance, 
causality and time-reversal invariance of electromagnetism as well as the idea 
that the Planck- Wheeler length (10 -33 cm) is the shortest scale allowable in 
any physical theory. The very generality of its assumptions guarantee the 
applicability of the scheme to other "gauge" interactions such as the strong 
forces. Besides, it introduces a useful simplifying assumption; namely, that 
the gravitational sector is unaffected by the scalar field introduced to vary the 
coupling constant . However, it is interesting to explore first this simplified 
model, before a similar exploration of more general theories. 

In fact, a grand unified description of electromagnetism and QCD will, 
presumably, predict linked time variation of low energy fundamental con- 
stants Nevertheless, one can anticipate that the QCD vacuum 
contribution, as we shall see in section 3, will be several orders of magnitud 
larger than the contribution coming from a em variability Hl3 |, so that we 



can assume, in the presence of a variable unifying coupling constant, that 
the strong forces separate cleanly from other interactions and so study them 
alone. To date we know of no "gauge-principled" analysis for the variability 
of the strong coupling constant and the object of this letter is to provide just 
such a study. In particular we will apply Bekenstein scheme to QCD since 
its assumptions are still valid for the strong interactions and find, contrary 
to the case of electromagnetism, that it is the vacuum, when compared to 
matter, which plays the dominant role as a source of variability. 
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2 Analysis 



Our starting point is the QCD Lagrangian with a varying coupling "constant" 

Lqcd = L e + L g + L m 

= L e - \TriGTGp,) + £ ^ f \i^d, -M f + g t{x)A^)^f\l) 
1 f 



where = G^t a , A* = A%t a , [t a ,t b ] = if abc t c , M f is the /-flavour 
quark mass and where, follwing ||, we introduced a classical scalar gauge- 
invariant and dimensionless field e(x). The varying coupling constant is given 
by g(x) = goe(x) where g Q is a constant and we require our theory governing 
e to be scale invariant. 



In order that the action be invariant under ip — > ip' = Utp (U = e 



-it a e a {x) 



) 



y._ ' .. {d,jU)U 1 while the gluon tensor 



we find that eA^ -> eA'^ = UeA^U' 1 - 

field, transforming like G — > G' = UGU~ l , should be given by 

1 

e 

Similar to the electromagnetic case, it is the time reversal invariance which 
excludes the G*G term from the free gluon action while, concerning the 
dynamics of e, the same arguments in || apply so we take 



G% = -[d,(eA a u )-d u (eAl)+g e 2 r bc A b ^A c u } (2) 



(3) 



2 I 2 e 2 

where we merely require the scale length I to be no shorter than the Plank- 
Wheeler length L P = y^f. 

Writing the Euler-Lagrange equations for the total action S = J LQCDy/—gd 4 x, 
first with respect to A^, we find 



" a 



• l-i 



gof^G^A^ + ^gajt^ 
f 







(4) 



then with respect to e, we get 



r 2 



□ In e = — 

he 



{G^Al), v + g -r bc G^A b ll A 



f 



(5) 
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Substituting (J|) in (|) we obtain 



□ lne 



he 



de * 2 a >' 



[IV 



(6) 



This equation is analogous to the electromagnetic case but while, in elec- 
tromagnetism, only matter acted as a source for both terms, by contrast, in 
QCD we can drop the first term working in the chiral limit Mf rs and, 
more importantly, we should include the vacuum contribution to the energy 
density in the second term. Approximating the mass in the universe by free 
nucleons we can write 



12 

□ In e = — 

he 



5<o|GTG>l°>- £ 



nucleons 



(N 


Tr(G^G^) 


N) 


(N 


N) 



2hc 



(0 I G^G a | 0) + Pm Af 



(7) 



where p m is the matter density of the universe and where the matrix element 
A)P is defined for the twist-2 operator by 



(N | Tr(G^ u G^ 2U ) - traces \ N) 



traces) 



and has the physical meaning of the part of nucleon momentum carried by 
gluons. 

Assuming homogeneity and isotropy for an expanding universe, and so 
considering only temporal variations for as, we get 



a 3 (t)l 2 c 2 
2hc ° 



(G 2 ) + Af Pri 



(9) 



where a(t) is the expansion scale factor. Let us consider models where the 
scale factor of the Universe behaves as in a flat Robertson- Walker space-time. 
Thus, for models with cosmo logical constant, we have 



a(t) = a(t ) 





r 3 i i 


\tt A J 


sinh -Q^H t 



(10) 



where f2 mj A is the cosmological density parameter corresponding to the mass 
and the cosmological constant respectively and Q m +Q\ = 1 for flat universes. 
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Models without cosmological constant are obtained in the limit Oa 
where the scale factor behaves as: 



a(t) = a(t Q ) 



0, 



(11) 



3 Results and Conclusion 



The matrix element can not be computed perturbatively and experi- 



ments give it the value of 0.48 14] while for the other non-perturbative quan- 
tity: the gluon condensate, we can estimate it by QCD sum rules method 



(0 | —G 2 | 0) ~ (0.012 ± 0.004) GeV 4 

7T 

where the operator G 2 is renormalized at the "natural" scale 1 GeV cor- 
responding to the matching condition of the sum rules to avoid the ap- 
pearance of large radiative corrections. In order to evaluate (G 2 ) we take 
h-QCD = 125 ± 25 MeV consistent with the range of values used in QCD sum 
rules flTj| implying, to leading log, ag(l GeV) = 0.336 ± 0.0323 and so we 
find 

(G 2 ) ~ (0.112 ±0.048) GeV 4 



Now since, in our model, p r , 



3H%c 2 



(^) 3 with Pc 



3ffgc 2 



10 



-47 



GeV 4 



8ttG \ a(t) J u rc 8nG 

we find, here, a strikingly exotic predominance of the QCD vacuum over 
matter and we can neglect the mass term altogether in equation @. By 
integration we find then for models with cosmological constant: 



X 



(t 



p c 16tt a 3 (t) 
t c ) ^ smh(3y/Q\H t) 



A 



siah.(3y/n^H t c 



;i2) 



where t c is an unknown free parameter. In the limit of zero cosmological 
costant, the variation of e behaves as: 



,2 (G 2 ) H 2 




(13) 
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In order to evaluate equations ( |I2"D and ( pT3| ) for today we need laboratory 
bounds on the variation of the strong coupling constant. Here, we can use a 
large number of data from various high energy processes ordered by increasing 
energy scale of the measurement as follows: r decay, GLS sum rule, QQ 
lattice, deep inelastic scattering, i?(e + ,e~), Pt(w), e + e~ event shape and Z 
width, giving in all a weighted average to Aq C d equal to 195 ± 65 MeV in 



the year 1994 |7| and 208 ± 25 MeV in the year 1999 |LJ. This allows us 
to take, up to leading log terms, the laboratory bound | - \today= 2 I ~ l< 
4.10~ 2 yr~ 1 - 

If we assume, plausibly, t c of the order of to ~ 10 10 yr ~ Hq 1 and take 
in the interval [0.25, 0.75] with Q m = 1 — Qa for models with cosmological 
constant (equation [L(|) or use equation (|TTD for models without cosmological 
constant, we find, pushing / down to near its minimum allowable value j- ~ 
1, a constraint on | t — t c \< I0~ 25 yr which is highly strange barring a "fine 
tuning" situation. 

For the purpose of refining the analysis, let us substitute equation (|i~0|) 
into (0) and integrate with e(to) = 1 to get the following expression for the 
variation of a s for models with cosmological constant 



= ^— 7^— [x coth x — xq coth x — x c coth x + x c coth xq 

a s \LpJ p c 127ri2 A 

+ sinh x c cosh x c coth x — sinh x c cosh x c coth Xo] (14) 
where x = ^y/TlAH t and xq(x c ) is x evaluated at t (t c ). In the limit of zero 



2 

c cosmological constant, we find using equation (|1 



A«5 (±Y(£^Hl(t_tl tl_tl\ 
a s \Lp) p c 8tt \ 2 2 t t J { ' 

Now, we can use astronomical and geophysical data giving bounds on the 
variation of as ranging over longer periods of time. In fact, the authors of 
[0] have derived a relation between the shift in the half-life of three long lived 
13 decayers ( 187 Re, 40 K and 87 Rb), measured in laboratory or by comparison 
with the age of meteorites, and a possible temporal variation of the funda- 
mental constants a em , Aqcd and Gf- In this work we attribute the change 
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uniquely to Aqcd and so we get a bound for the variation of as at the age 
of the meteorites compared to its value now = (0 ± 2.1 x 10~ 4 ). 

One of the most stringent limits on the time variation of fundamental con- 
stants follows from an analysis of isotope ratios of 149 Sm/ 147 Sm in the natural 
uranium fission reactor that took place 1.8 x 10 9 yr ago at the present day site 
of the Oklo mine in Gabon, Africa [[19], 20| . Sisterna and Vucetich || have 



shown that the information about the time variation of the strong interac- 
tion parameter Aqcd is very small when compared to the other fundamental 
constants. This is because in the chiral limit, the time variation of any strong 
energy difference takes the form ^Jj = 4; and the limits on the shift of the 
neutron capture resonance are not very accurate. (See also the discussion in 
[^IJ.Thus, we will not include the Oklo limits in our set of data. 

On the other hand, quasar absorption systems present ideal laboratories 
to test the temporal variation of the fundamental constants. The continuum 
spectrum of a quasar was formed at an epoch corresponding to the redshift 
z of the main emission details with the relation X Q b s = Xi a b (1 + z). Knowing 
that the ratio of frequencies of the hyperfine 21 cm absorption transition 
of neutral hydrogen to an optical resonance transition is proportional to 
x = o? em 9 v ^ where g p is the proton g factor and m p is its mass, we can 
translate a change in x into a difference between the measured redshifts of 
the 21 cm and the optical absorption as follows: 

Ax z opt — Z 2 i , > 

— = w ( 6) 

Thus, combining the measurements of optical and radio redshifts, one can 



obtain bounds on x: = (0.7 ± 1.1) x 1(T 5 at z = 1.776 pi 
(0 ± 1.2) x 10~ 4 at z = 0.69 |!| and ^ = (0 ± 2.8) x 10" 4 at z = 0.52 |3 



Ax 

x 



Evidence for the time variation of the fine structure constant has been 



claimed by Webb et al 0] and further results posted recently [gq, 26| . How 



ever, this detection does not suggest a power law fit. Moreover, it is found in 



27] that Webb's results are in conflict with the equivalence principle when 



examined within Bekenstein model for a em variability which is in agreement 
with the conclusion of [13| examining all available experimental bounds in 
the context of typical theories predicting time variation of the fundamental 
constants. For the purposes of this paper, the bounds on ^-variations will be 
assumed to come solely from the change of m p proportional to Aqcd and, 
thus, infer bounds on as-variations. Moreover, observations of molecular 
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hydrogen in quasar absorption systems can be used to set bounds on the 
evolution of fi = — throughout cosmological time scales: — < 2 x 10~ 4 at 



z = 2.811 i28| , and this, in turn, would imply a bound on as-variation under 
a similar assumption to that for the other quasar data. Bounds on the mass 
ratio have also been examined (p9l [30| using comparison of 21-cm hyperfine 



hydrogen and molecular rotatinal absorption lines. However, it is shown that 



C) 



this comparison actually constraints a em g p pTT 

Taking the data described above, with more "reasonable" choices for t, 
we have performed a statistical analysis using MINUIT with H Q = 6.64 x 
10~ n yr _1 and obtained the results in the table [l] for the free parameter in the 
model y- with 99 % of confidence level. Furthermore, we could check that, 
as long as t c is far from the unnatural "fine tuning" situation, we always get 
/ shorter than Lp by too many orders of magnitude in clear conflict with the 
postulates adopted in the framework. While integrating, we have assumed 
the value of the gluon condensate < G 2 > constant whereas it might vary 
over cosmological time scales. However, we do not expect such variations to 
be large enough to change the conclusion above. Even though length scales 
shorter than the Planck- Wheeler length Lp might enter physics in the context 
of "new" theories, namely superstrings, it is extremely unlikely that our tiny 
I could be accommodated so as to recover the axioms of the Bekenstein 
model. Since the assumptions of this model are reasonable, comparison with 
experiments, excludes any direct generalization of the model to QCD and 
consequently rules out a s variability in accordance with the strong principle 
of equivalence. Although it is possible that more general schemes, involving 
more realistic assumptions on the structure of the e field or its interactions 
with gravitation, may resist comparison with the experiment, the present 
study shows that it will be difficult to build them from first principles without 
some fine tunning mechanism. 
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Table 1: Constraints on 4- using the full data set of bounds on for flat 
models (fl\ + fl m = 1) with different cosmological constant. t c = jto 





7 


j_ 




0.75 


0.3 


(0±5) x 10' 


-25 
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(1 ±4) x 10" 


-25 
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(0±6) x 10' 
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(2 ± 7) x 10- 
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(0 ± 7) x 10- 
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(4 ±8) x 10- 


-25 
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(0 ± 7) x 10- 


-25 
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